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ABSTRACT 

The general model for incoherent synchrotron radiation has long been known, with 
the first theory being published by Westfold in 1959 and continued by Westfold and 
Legg in 1968. When this model was first developed it was applied to radiation from 
Jupiter, with a magnetic field of ~ 1 G. Pulsars have a magnetic field of ss 10^^ G. 
The Westfold and Legg model predict a circular polarization which is proportional to 
the square root of the magnetic field, and consequently predicts greater than 100 per 
cent circular polarization at high magnetic fields. Here a new model is derived based 
upon a more detailed analysis of the pitch angle distribution. This model is concerned 
with the frequency range /bo/7 << / ^ /boj noting that /b„ = 2.7 x Iff B, which 
for a relatively high magnetic field 10® — 10® Gauss) leaves emission in the optical 
range. This is much lower than the expected frequency peak for a mono-energetic 
particle of 0.29 ^®^^ 7^. We predict the circular polarization peaks around lO^G in 
the optical regime with the radiation almost 15 per cent circularly polarized. The linear 
polarization changes from about 60 to 80 per cent in the same regime. We examine 
implications of this for pulsar studies. 
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1 INTRODUCTION 


Pulsar emission in the optical regime is generally accepted to be incoherent synchrotron radiation and consequently it should be 
polarized. To date most a t tention has been o n linear polarization, in part due to instrumental limitations of most polarimeters 


Slowikowska et al 


Collins et. a' 


( 2009 1^ Smith et al ( 1988ll . Optical instrumentation, such as the Galway Astronomical Stokes Polarimeter 
,]| ( 201,1 1. are now in a position to measure all of the Stokes parameters from pulsars on time-scales from 


(GASP) 

milliseconds to hours. Hence the requirement for a fully self-consistent model for synchrotron radiation in a high magnetic 
fields. 

The original mo del for synchrot r on emission was pnbl is hed in 1959 WestfokJ Jlj lSgf), with other authors coming to 
the same conclusions Le Roux 1 1961 1: Ginzburg et al. ( 1968l l: Ginzburg fc Svro vatskiil (|l96E ), al b eit ge nerally with slightly 
different derivation methods. T his model was then f urther developed in 1968 Legg fc Westfold ( 1968h . and corrections to 
the model were applied in 1974 iGleeson et al.l lll974h . Thes e corre c tions do not significantly alter the circular polarization 
calculations. An error in the derivation used was found 1986 Singall 1 1986h . but this did not change the model predictions. 

The motivation behind the Westfold and Legg (hereafter WL) model of synchrotron radiation was to study the emission 
from Jupiter, with a magnetic field of approximately 1 G. As such, the behaviour of the model was never tested at high 
magnetic fields. As pulsars have extremely high magnetic fields (~ lO^^G), it is important to test the model in high magnetic 
fields before applying it to pulsar emission. When the incoherent synchrotron emission is calculated at high magnetic fields 
(~ 10® G), the WL formulation predicts a circular polarization greater than 100 per cent. As this is in clear contradiction of 
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reality, a new model for the incoherent synchrotron emission is required. The most likely reason was in the expansion of the 

particle pitch angle distribution. In this paper, this is expanded to the next order of magnitude. _ _ 

The second section of this paper states the predictions of the previous model Legs fc Westfold ( 1968h . and gives the 
problems inherent in that model. The third section goes through the expansion of the particle pitch angle distribution and 
the differences that this makes to the intensity, linear polarization and circular polarization. Finally, the overall effects of each 
of the different parameters on the intensity are investigated, and some of the possible implications are discussed, particularly 
in relation to pulsar emission. 


2 THE WL MODEL OF SYNCHROTRON EMISSION 


The WL model that is of interest here gives the Stokes parameters for a power-law distribution of electrons with 
N{E) = E-^ Ei< E <E2 

and N{E) = 0 otherwise as 
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where / is the frequency, p is the permeability constant, p is the power-law index, e is the charge of the particle, 9 is the 
particle pitch angle, fso = eB{2nmc)~^ is the fundamental gyro-frequency, $(0) is the pitch angle distribution function, and 
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The circular polarization (defined as VI is then given by 
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As the circular polarization is proportional to the root of the magnetic field, it is clear that at some point the degree of circular 
polarization will exceed 1. This is clearly unrealistic. However, all models will only have a certain range of validity, and if the 
magnetic field at which the circular polarization occurs at extremely high magneitic fields, then the model can still be used 
for smaller magnetic fields. The polarization was found to be greater than 100 per cent at approximately 10® to 10^ gauss 
(Fig. [ 1 ]) and above, below the surface magnetic field strength of pulsars, but above planetary magnetic fields. 


3 EXPANDING THE WL MODEL 

As the model fails at high magnetic fields, a new model is needed to describe the polarization of synchrotron radiation in those 
fields. When the derivation used in the WL model is looked at, the main assumptions are 1) that the velocity of the particle 
is close to the speed of light 2) the expansion of the velocity in the frame chosen 3) for a power-law index, the expansion of 
the distribution of electrons. 

In this work we will examine the results when the distribution of electrons is expanded to a higher order. 


3.1 Electric Field 


As synchrotron radiation comes from a source moving in a cyclic fashion, emission will consist of harmonics of the fundamental 
gyro-frequency /sp. Following the same formulation as WL, the emission from each harmonic can be shown to be 

n X ((n X X 
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% Circular Polarisation 

Frequency = 5.2 x Hz 
p = 1.42 



Figure 1. The WL model predicts that the percentage circular polarization will increase linearly with the magnetic field, regardless of 
the power-law index of the electrons used. At some point this model fails and predicts clearly non-physical results. This is dependant on 
the pitch angle, the frequency, and the power-law index, but in the area of interest the WL model fails long before the predicted surface 
magnetic field of a pulsar 10^^ Gauss). 


where the expression in brackets is evaluated at the retarded time 
R{t') 


t =t- 


R(t^) = r — r{t') 


r — n • r 


where b = /3'sinasin(a — 9) Le Roux ( 196lh : Legg Sz Westfold 

over t', and simplifying, gives 


E^ = 


/ice 

Svr^r 


(x)' *^exp (in {^) I n X (n X /3) exp [m (t' - 


(9) 

. Changing the integration to an integration 

dt' (10) 


3.2 Coordinate System 

In order to solve equation cni, it is necessary to choose a system of coordinates. In this case, the system will be constructed 
as follows: the particle is spiralling around a magnetic field at an angnlar freqnency of ujb = qB{'ymc)~^, where q is the 
charge, B is the strength of the magnetic field and 7 is the Lorentz factor of the particle. The particle maintains a constant 
pitch angle of a with respect to the the magnetic field direction. At any particnlar time the orbit has a radius of curvature 
of a. Now, let the x-y plane be the instantaneous plane of the orbit of the particle. Now, take the origin of the a;-axis to be 
the point where the velocity vector and the observer are in the x-z plane, and let the y coordinate be in the direction of the 
radial vector a, with the x coordinate being defined as perpendicular to the y and 2 coordinates. 

Now, define a new set of coordinates (n, ey, ej_) such that the origin is at the same point as the {x,y,z) coordinate system 
origin, n is pointing towards the observer, e± is pointing along y, and ey =nXe±. This then gives a natural coordinate system 
























4 D. de Burca & A. Shearer 



Figure 2. The geometry used in order to calculate the synchrotron emission. Firstly, define the x-y plane as the instantaneous plane of 
orbit of the particle. Then define the origin as the point at which the velocity v and the vector to the observer n are both in the x-z 
plane. Define ey to be along the y axis, and ey as nXey. This gives a natural frame of reference for the polarization of the emission. 


in which to consider the polarization of the emission, as ey is perpendicular to the magnetic field and ey is parallel to the 
magnetic field direction, as seen in projection by an observer, as can be seen in Fig. O Here a stands for the radius of 
curvature of the particle. _ _ 


Finally, in this coordinate system the velocity and acceleration are iLongairi ll201ll ~) 
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This then gives the electric field (splitting it into its component parts parallel and perpendicular to the projection of the 
magnetic field), and dropping the subscript, 
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3.3 Emission-polarization tensor 

The emission polarization tensor is defined as 
27rr^ / Ea_E*^ ExEy \ 

^ n \ E^EI E||Ei| j 

This is equivalent to getting the Stokes parameters for each harmonic, as 
I = pil + P22 
Q = Pll — p22 

U = Pl2 + P 2 I 

V = - (pi 2 — P 21 ) 

I 


(19) 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


3.4 Airy Functions 

It is possible to convert the electric field exponential into Bessel functions (see supplementary materials). This gives 
, 4 
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This then gives the polarization tensor for a particular harmonic of the emission. 


(24) 

(25) 

(26) 


3.5 Converting to Frequency Domain 

For large-order harmonics, the radiation becomes quasi-continuous Legs &: Westfold (196^) and it is possible to convert the 
polarization tensor for a single harmonic to the frequency polarization tensor using 
Pub 


PS = 


/b 


(27) 


where /s = /sq/Ti b = /?'sinasin (a — 6^), pn is the polarization tensor for a single harmonic and p/ is the polarization 
emission tensor at a particular frequency, and 

/ = (28) 

This gives the polarization tensor at a particular frequency. It is convenient to convert from the frequency into a dimen¬ 
sionless parameter x such that 
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3.6 Power-law Polarisation-emission tensor 


When there is a power-law of particles, the polarization-emission tensor for that population of particles is 
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In order to solve this it is possible to represent a as a-\- 6. Then the solid angle is represented as df2(n) as 27rsinad0. This 
can be written as a = a' + 6 and substituted into equation dSl. To third order, the particle pitch angle distribution can be 
written as 

4>{a + 6) sin [a + 61) = fia) + g{a)9 + h{a)9^ (35) 

where 

f{a') =(fi{a)sma' (36) 

g{a) =<j)'(a) sin a' + <j){a) cos a (37) 

h{a) =(j>"{a) sin a + 2(f)'(a ) cos a — (j>{a) sin a' (38) 

and (f){a) is the pitch angle distribution of the particles itself. This gives the polarization-emission tensor as (writing a' as a, 
see supplementary materials) 
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and the other parameters are as before. 


4 DISCUSSION 


This formulation predicts that the cirular polarisation will not exceed one hundred per cent. For a particle power law index 
of 1.42, the circular polarization remains less than fifteen percent for all magnetic field values in the optical regime (Fig. 
o, while the linear polarization changes value at high magnetic fields (Fig. [4]) but remains at a steady value except in the 
intermediate range of magnetic fields (10® — 10^® G). One result is that the level of circular polarization is directly dependant 
on the magnetic field. As such, measurements of the circular polarization gives some potential constraints to the magnetic 
field strength of the emission volume. 

There is one major constraint to this emission model. In order to obtain emission at high magnetic fields, the h{a) term 
has to be positive. In general, for small values of 9, this can be obtained by considering particle pitch angle distributions with 
positive first derivatives. So, in this paper the particle pitch angle distribution used (f>{(y.) = sin{a)sin{amax)~^■ However, this 
particle pitch angle distribution does not in general agree with the particle pitch angle distributions that are predicted (e.g. 
guassian). Physically, this type of distribution would result from particles which lose more energy the closer they are to the 
magnetic held line. 

Another area of interest would be in pulsa r studies. The r e are c urrent l y a n u mber of d ifferen t theories about pulsar 
high-energy/optical emission Cheng et a] 1 1986ll : Takata et al 1 2O10h : Petri (2012); Du et al ( 2011 1: Harding ( 2013ll which 
agree on the process of emission, pair production creating a plasma which then emits using synchrotron radiation, but which 
disagree on the location of the pulsar emission zone. As the magnetic held strength can be correlated with the position in 
the pulsar magnetosphere, this provides a test to constrain the p ulsar emission location. One method to constrain the pulsar 
emission locations would be to use an inverse mapping approach. iMcDonald et al.l (1201111 . They considered emission from all 
parts of the magnetosphere, and compared that emission to optical observations. They found that the majority of the emission 
came from approximately 300 km from the pulsar surface, where the magnetic held strength is in the range of 10^ — 10® G. 
Future work would involve incorporating our model into the code and checking if the emission areas change signihcantly. 

The linear polarization predicted by our method is not dependant on the particle pitch angle except in intermediate 
magnetic helds, and in certain regimes is not dependant on the magnetic held strength. However, it is very sensitive to 
the particle power-law index. As such, it could be possible to constrain the particle power-law index from measurements of 
the linear polarization. A combined measure of linear and circular polarization is therefore an important diagnostic tool for 
determining the geometry of pulsar emission zones. 

The relationship between the observed power-law index , and the partic l e power-law index, is d ifferen t at high magnetic 
helds. This relationship has been accepted as p = 2a -|- 1 Long:aii ( 2011 1: iRvbicki fc Lightman ( 1979l l. however, at high 
magnetic helds, this relationship changes to p = 2 q — 1 in this formulation (Fig. [5]). 
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% Circular Polarisation 

Frequency = 5.2 x 10^“^ Hz 
p = 1.42 



Method 



Figure 3. The circular polarization for a power-law distribution of particles with a power law index of 1.42, at a frequency of 5.212 X 10^^ 
Hz. Here WL stands for Westfold and Legg, the original emission theory, and theta stands for the particle pitch angle. 


There are still a number of issues to be addressed. The frequencies of interest here are far from the maximal spectral 
frequency of a single shape. Therefore, it is worthwhile to note that there conld be errors introduced due to the integration 
going between 0 and oo rather than over a realistic energy range. Investigating the error is beyond the scope o f this work. To 


first o rder we can state that the polarisation measurements are correct, as calculating the formula for and .7„ iGleeson et al, 

( 197J) gives a ratio of 0.644 for the integration over 0 to oo , and a slowly varying ratio between 0.644-0.500 (as the magnetic 
field is increases from lO^-lO^*^) for the ratio from the exactly calculated values. Lifetime effect s also limit the effecti ve energy 
range - if the energy is below a minimum of 7 = fbo/if sitd 0) then there can be no radiation lOleeson et al.l 1 19741 ). At high 
magnetic fields this can be of the order of 7 = 10 ®. 

Measurements of the linear and circular polarization from pulsars with apparent magnitudes less than 25 is possible 
with instruments like GASP on 4-m class telescopes. Our predictions can therefore be tested on normal pulsars such as the 















8 D. de Burca & A. Shearer 


% Linear Polarisation 

Frequency = 5.2 x 10^"^ Hz 
p = 1.42 


80- 



1 100 10 ® 10 ® 10 ^° 10 ^^ 

Magnetic Field(Gauss) 


Figure 4. The linear polarization change with regard to the magnetic field for a particle power law index of 1.42 and frequency 5.212 X10^^ 
Hz and theta for the particle pitch angle.. As can be seen, the linear polarization is steady at low magnetic fields and at high magnetic 
fields, with the linear polarization changing smoothly between the two values in intermediate range of magnetic field values. 


Crab pulsar and on magnetars snch as 4U0142+61. We also develop onr inverse mapping approach [McDonald et al 
to inclnde circnlar polarization. 


(l201lll 


We note that for the Crab pulsar the maximum linear polarisation is fifteen per cent Slowikowska et al. 1 20091 ') . whereas 


we predict higher values (> 60 per cent), consistent with the WL formulation. This discrepancy requires further investigation 
and could be due to either the impact of different pitch angle distribution or a more astrophysical explanation. Future work 
should clarify this. 
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Intensity Vs Magnetic Field 


Frequency = 5.2 x 10^'* Hz 
p = 1.42 



theta 

— 0.01 
.... 0.1 
0.2 
- - 0.5 


Method 

^WL 
— This Work 


Figure 5. A comparison of the intensities predicted by both models. The frequency is 5.212 X 10^^ Hz, the pitch angle distribution is 
il?(Q:) = sin(«), the particle power-law index is 1.42, and theta stands for the pitch angle. As can be seen, at approximately 10^ — 10® G, 
the slope of the intensity changes from a = 0.5(p — 1) to 0.5(p + 1). 
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